Vortex matter and generalizations of dipolar superfluidity concept in layered systems 
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In the first part of tiiis fetter we discuss efectrodynamics of an excitonic condensate in a biiayer. 
We sliow that under certain conditions the system has a dominant energy scaie and is described 
by the effective efectrodynamics with "pfanar magnetic charges" . In the second part of the paper 
we point out that a vortex fiquid state in bifayer superconductors afso possesses dipofar superfluid 
modes and estabfish equivafence mapping between this state and a dipofar excitonic condensate. We 
point out that a vortex liquid state in an N-fayer superconductor possesses muftipfe topofogicaffy 
coupled dipolar superfluid modes and therefore represents a generalization of the dipofar superflu- 
idity concept. 



The progress in semiconductor technology has made it 
feasible to produce bilayers where the interparticle dis- 
tance is larger than the separation of the layers and 
which are expected to have interlayer excitonic states 
P, S, H, i, il (as shown in Fig. 1) with a significant 
life time. Some signatures of Bose-Einstein condensate 
of such excitons were reported_[6||. A renewed theoretical 
interest in these systems [H, [3, 01 is focused on identi- 
fication of possible unique properties of such condensates 
which can set them apart from other families of quan- 
tum fluids. In [4] it was discussed in great detail that the 
spatial separation of the positive and negative charges, 
in an exciton in a biiayer, makes the phase of the con- 
densate transform as V0 \/9 — eA(r+) + eA_(r_). 
The dipolar excitonic condensate (DEC) therefore fea- 
tures a coupling to a difference of vector potential val- 
ues at positions in different layers r_(- and r_ P, [3, 13|- 
In the case of a small layer separation, d, one finds 
A(r+) — A_(r_) « ddzA{r). Therefore a static in-plane 
magnetic field H*^^* produces excitonic currents and the 
system is described by the following free energy density 



I (ye + ed[RW*]y 



(1) 



where the operator R rotates a vector 90 degrees coun- 
terclockwise: R{a,b,c) — {—b,a,c), p is the phase stiff- 
ness, e is the electric charge and d is the separation of 
the layers. A particularly interesting aspect of this ob- 
servation is that although an exciton is, by definition, 
an electrically neutral object, it has the dipolar coupling 
to a gauge field and its effective low-energy model has a 
symmetry different from the so-called "global" U{1) sym- 
metry which one finds in ordinary neutral supcrfluids. It 
is also different from the "gauged" U{1) symmetry which 
one finds in superconductors. 



To describe topological defects we should include dynam- 
ics of the gauge field in the model ([3]). Consider a sys- 
tem of two thin parallel layers with positive and negative 



FIG. 1: (color online) A schematic picture of DEC. An ex- 
citon forms as a result of the pairing of an efectron in an 
eiectron-rich iayer and a hofe in a hoie-rich iayer. 



charges bound in pairs (as shown in Fig. 1). Here, 

V X B = J = e[(5(z_) - 5{z+)] x 

|1(^V^* - ^*VV^) - e\yj\'[A{r+) - A(r_)]} (2) 

where J is the electric current, tp is the DEC order pa- 
rameter, _) are the z-axis positions of the upper 
and lower planes and r(-i-,-) — (x, ?/, Zj-^ ^-j). Topolog- 
ical defects in this system correspond to the situation 
when the phase, 6, of the order parameter, ip, changes 
by 27m. A configuration of accompanying magnetic field 
should be determined by minimization of the energy tak- 
ing into account (i) the kinetic energy of currents in 
two planes, (ii) the potential energy of the corresponding 
Ginzburg-Landau functional, (iii) the energy of the three 
dimensional magnetic field configuration. This problem 
is non-local and the field-inducing current itself depends 
on the gauge field, i.e., should be determined in a self- 
consistent way. However we show here that there is a 
regime when the system is accurately described by an un- 
usual, on the other hand tractable effective model. That 
is, let us consider the situation where the separation of 
the layers d is small compared to the system size and 
A(r+) — A_(r_) « ddzA{r). Then in the hydrodynamic 
limit the effective model is 



V0 + ed[RB,n{v)] 



-B'(r) 

2 ^ ' 



(3) 



Where Bi„ is the field in the dipolar layer. Consider a 
vortex with A9 — 2tt. Then V0 = ieg, which produces a 
logarithmic divergence of the energy in a neutral system. 
However eq. ([3]) suggests that for a given phase winding, 
the system can minimize its energy by generating a cer- 
tain configuration of Bi„. Note that the second term in 
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([3]) depends quadratically on B and does not allow a con- 
figuration of Bin which would completely compensate the 
phase gradient in the first term. A configuration of the 
interplane field, which would partially compensates the 
divergence caused by V6', should satisfy the condition: 
ec?(i?Bi„) oc (r 3> Tcore)- This implies the following 
self-induced interplane field: 

1 r 

B,„ = a— — [a < 0,r = (x,y,0)] (4) 
ed 

At a first glance, eq. ([4|) appears to violate the con- 
dition that the magnetic field should be divergenceless. 
This problem is resolved by including in the picture the 
"out-of-plane" magnetic fields. These, as schematically 
shown in Fig. 2, can restore the V • B = constraint 
in three dimensional physical space while permitting the 
in-plane field to have the form required by ([4]) with a nat- 
ural cutoff close to the vortex core. From the V • B = 
condition it follows that both the interlayer Bi„ field and 
the out-of-plane field carry the same flux. However for a 
given magnetic flux, the out-of-plane field has the free- 
dom to spread in the positive and negative z-axis direc- 
tions. Since the magnetic flux is J B • dS, while the mag- 
netic field energy is (1/2) J B^, the out-of plane field will 
have a finite value of the integral (1/2) /^.^^j^ z+]'^'^ 
over entire three dimensional space excluding the bilayer 
space. On the other hand the magnetic field inside the 
bilayer behaves as |Bi„| oc 1/r and thus has logarith- 
mically divergent energy (1/2) /^^gj^ Therefore 
from the condition V-B = and geometry of the problem 
it follows that the energy of out-of-plane magnetic fields 
is negligible compared to the energy of magnetic field in 
the dipolar bilayer for a sample with d small compared 
to the system size. 

Thus we have identified the regime where the dynam- 
ics of the magnetic field is dominated by its most ener- 
getically costly (weakly divergent) interlayer part, Bi„, 
which leads to an interesting two-dimensional effective 
model where the magnetic field Bi„ is not subject to the 
constraint that V • Bi„ — 0: 



Side view 




FIG. 2: (color online) A schematic picture of a magnetic field 
configuration in DEC with vortices. Side view: the out-of- 
plane magnetic field B carries the same flux as the interplane 
field 'Bin, however it has the freedom to minimize energy by 
spreading above and below the top and bottom layers. Top 
view: the in-plane magnetic field Bi„ configuration for vor- 
tices with phase windings A9 = ±2n 



5d[RBin{x,y)] 



-Bi(x,y)(5) 



Let us consider a vortex with a phase winding A6 — 27rn 
in the model ([5]). The coefficient, amm, for the ansatz 
Q which minimizes the spatially integrated free energy 
density (O for vortex with A9 = 2nn is a^j^^'^" = 
—np[p+ (e^d)~^]~^. Thus the vortex has the following 
configuration of the magnetic field 



B 



Ae=2TTn 



P 



1 r 
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(6) 



where r — {x,y,0) (shown on Fig. 2). Therefore, 
these defects emit a quantized radial magnetic flux $ = 
B(r)|27rrd. The quantization condition is: $^^=2?™ _ 

e p+{e-^d)- 



^ s_i n. Thus a "dipolar flux quantum" is 



$2 



2tt pe^d 
e pe^d + 1 



(7) 



where $o = 27r/e is the standard magnetic flux quan- 
tum. In a way these topological excitations play a role of 
positive and negative "magnetic charges" in the effective 
planar electrodynamics of the model ([5|) . 

The energy of such a vortex placed in the center of a 
circle-shaped system of a radius R is: 



E 



■npn 



:ln 



R 



1 + pe'^d Tc 



(8) 



where r^ore is the vortex core size. Even though this 
energy is logarithmically divergent, these vortices can in 
principle be induced in a DEC by an external in-plane 
magnetic field because their magnetic field also has diver- 
gent energy. This provides a possibility to obtain a nega- 
tive contribution to the Gibbs energy from B ■ H. For 
example, in principle, an experimentally realizable con- 
figuration of induction elements can produce the follow- 
ing external in-plane field H w [r = {x, y, 0); r > 
ro] where tq is some cutoff length. Then, for a vortex 
with A9 = 27r, the integral B • H provides a negative 
logarithmically divergent contribution to the Gibbs en- 
ergy G — /^(£ — B • H) . There is a critical value of 
the coefficient, 7 at which the creation of a vortex be- 
comes energetically favorable (below we set tq w Tcore 
without loss of generality). In the regime pe^d <C 1 we 
obtain 7c « In general case the energetic favora- 
bility of a vortex state depends on the integral B • H. 
An external field which is significantly stronger than that 
corresponding to 7c favors a higher phase winding num- 
ber of the induced vortex structure, however on the other 
hand, the vortex energy depends quadratically on n [see 
eq. ([5])]. Therefore, a stronger field should normally pro- 
duce a state with multiple one-dipolar-quantum vortices. 
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In a planar C/(l)-symnietric system, thermal fluctua- 
tions can excite finite-energy pairs of vortices with oppo- 
site windings. In the model ^ the interaction between 
a vortex with A9 = 27r located at ri and an antivor- 
tex with ^ —2tt located at r2 originates from two 
sources. The usual current-current interaction produces 
the attractive force 



Fj = -27rp 1 



pe^d 



(r2 - ri 



pe^d +1 J |r2 - ri| 



the other contribution comes from the Bf„ term: 

pe'^d (r2-ri) 



-27rp 



[pe^d + 1]^ |r2 — ri 



(9) 



(10) 



A system of these vortices and antivortices can be 
mapped onto a Coulomb gas which at a temperature Tkt 
undergoes a KT transition: 



2 l + p(TKT)e^d 



(11) 



Therefore, the temperature of the condensation transi- 
tion in DEC will be suppressed compared to the value of 
the condensation temperature in a neutral system with 
similar density. While the suppression for realistic DECs 
is tiny, there is an interesting aspect in it because in 
contrast to the superfluid density jump p/Tkt = 2/7r in 
regular superfluids here one should define a "generalized 
superfluid density jump" which depends on the electric 
charge and on a nonuniversal paramter: the layer sepa- 
ration. 

As is well known, one of the ways to detect a KT 
transition/crossover in superconducting films is associ- 
ated with a peculiar reaction to an applied current Q. 
In a superconducting film an external current results in 
Lorentz forces acting on a vortex and antivortex in op- 
posite ways, causing a pairbreaking effect. Free vortices 
create dissipation which is manifested in IV character- 
istics 0. In DEC, vortices have magnetic field which 
can be viewed as that of "magnetic charges" and, cor- 
respondingly, they are sensitive to an applied external 
uniform in-plane magnetic field. Such a field, at finite 
temperature, should create a KT-specific modification of 
the zero-temperature response discussed in Q. 

We note that because of small carrier density and layer 
separation in the presently available semiconductor bilay- 
ers a dipolar vortex would carry only about 10^'' fiux 
quanta, which, though may be resolved with a modern 
SQUID, makes observation of these effects difficult. This 
raises the question if there could be strong-coupling dipo- 
lar superfluids in principle. Below we show that the con- 
cept to dipolar superfluidity arises in a system principally 
different from DEC, without interlayer pairing problem 
(which limits dipolar coupling strength in DEC). This 
provides a possibility to have larger ffux of a dipolar vor- 
tex and more pronounced phenomena associated with it. 



Moreover there the concept of dipolar superfluidity al- 
lows generalization. 

Consider a layered superconductor (LSC), i.e. a multi- 
ple superconducting layers separated by insulating layers 
(to effectively eliminate interlayer Josephson coupling) so 
that the layers are only coupled by the gauge field. This 
system has been extensively studied in the past 0]. In 
the hydrodynamic limit its free energy density is 

N ^ 

F = ^^S{zi) (v -ieA{x,y,Zi)j'ipi{x,y,Zi) 



(V X A{x,y,z)y 



(12) 



where 'ij;i{x,y,Zi) = |'(/'j(a;, y, Zj)|e'^*(^'^'''*\ 2;^ + = z^+i. 
This model indeed does not feature dipolar superfluidity 
of the DEC type. Here, the main distinction is the reac- 
tion to the external fleld which is screened because of the 
Meissner effect (the effective screening length in a planar 
superconductor is A^/ [layer thickness]). However there 
are situations where the superconductivity can be elim- 
inated in this system by topological defects. That is, if 
to apply an external fleld along ^-direction one can pro- 
duce a lattice of vortex lines with phase winding in each 
layer (A6li(zi) = 2tt, ...,Mn{zn) = 27r) d. Enclosed 
magnetic flux gives such a vortex line a flnite tension. At 
elevated temperatures the lattice of these vortices melts 
but importantly normally there is a range of parameters 
(which depends on the strength of the applied fleld and 
temperature) where the vortex lines forming a liquid re- 
tain tension [1, 0] . If the vortex lattice is not pinned or 
if one has a tensionfuU vortex liquid the charge trans- 
fer in superconducting layers is dissipative. However in 
such situations a system nonetheless retains broken sym- 
metries associated with the phase differences between 
the order parameters, and correspondingly dissipation- 
less countercurrents [lO| , in this particular case a broken 
symmetry is retained in the phase difference between the 
layers. Physically the situation which occurs is the fol- 
lowing: consider the TV = 2, |i/;i| = |?/'2| = IV"! case. Cur- 
rents in individual layers move unpinned vortices which 
produce dissipation. However as long as a vortex line 
threading the system has a flnite tension, equal counter- 
currents in different layers deform but not move a vortex 
line and therefore do not create dissipation. For small 
layer separation d, the dissipationless counterflows can 
approximately be described by extracting phase differ- 
ence terms [10|. Then the part of the model ([12]) which 
retains a broken symmery is: 



1 



V(6'i(a;,y,zi) - 6'2(x, za)) 



-e(A(a;,y,zi) - A{x,y,Z2)) 



(13) 



The vortices in this system are related to dipolar vortices 
in DEC. The simplest vortices with a topological charge 
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in the phase difference are (A^i = ±2tt, A02 = 0) and 
(A^i = 0, A62 = ±2tt). We denote them by a pair of inte- 
gers (±1, 0) and (0, ±1). For a vortex (±1, 0) the currents 
in two layers are: ji — elV'pV^i — e'^\'ip\'^A{x, y, Zi); j2 — 
— e^|'(/'|^A(x, ?/, Z2). Like in a DEC, in the case of a fi- 
nite d, the vortex features in-plane radial magnetic field 
[RBin] ~ dzA{r). These vortices can be induced by an 
external in-plane magnetic field, like in a DEC. However 
there are indeed also principal differences. LSC is a sys- 
tem with more degrees of freedom and the above consid- 
erations apply only to the state when superconductivity 
in individual layers is removed e.g. by a molten lattice 
of (1,1) vortices, or thermally excited (±1,±1) vortices. 
Note also that the molten lattice of (1,1) vortices does not 
automatically preclude a formation of an ordered state of 
(±1,0) and (0, ±1) vortices because a vortex (1,1) does 
not have a topological charge in the phase difference sec- 
tor V{9i — 62), and their disordered states can not ehm- 
inate corresponding phase stiffness. The density of (1,1) 
vortices and correspondingly the temperature of their lat- 
tice melting is controlled by the strength of the magnetic 
field in z-direction while the density of (±1, 0) and (0, ±1) 
vortices is controlled by the in-plane magnetic field. The 
system therefore possesses a control parameter which al- 
lows for the ordered structures of (±1, 0) and (0, ±1) vor- 
tices to coexist with a liquid of (1, 1) vortices. 

The dipolar superfluidity in LSC has a number of de- 
tectable physical consequences. Namely the LSC in the 
vortex liquid state should have a dipolar superfluid re- 
sponse, analogous to that discussed in great detail in Q 
. Also the system in the vortex liquid state should possess 
aspects of planar electrodynamics with effective magnetic 
charges discussed above in connection with DEC which 
may be observable in the temperature dependency of the 
dipolar response. The general case of N layers, especially 
with the variable interlayer distances and the condensate 
densities, has much richer structure than DEC because 
the dipolar superfluid modes are multiple and coupled. 
That is, for N layers one can have multiple combina- 
tions of counterflows in different layers which will not 
move a tcnsionfuU vortex line, therefore one has to con- 
sider all possible combinations of phase differences to de- 
scribe dipolar modes. In the limit of very small d the 
kinetic terms for countercurrents approximately can be 
expressed as: 



-e{A{x,y,Zi) - A{x,y,Zj)) . (14) 



tions will in some cases remain Q. 

In conclusion we have considered the possible physi- 
cal effects dictated by the symmetry and dynamics of a 
gauge field in dipolar condensates. For DEC, we started 
by constructing an effective model, based on symmetry 
and energy scales of the problem. In this framework we 
described topological defects which emit a nonuniversally 
quantized radial magnetic flux. We pointed out that be- 
cause of the existence of well separated energy scales in 
the regime, when interlayer distance is smaller than other 
length scales in the problem, the system possesses ef- 
fective planar electrodynamics with "magnetic charges" 
(planar analogue of magnetic monopoles) arising from 
topological excitations. Therefore at finite temperature 
an applied inplane magnetic field should have a vortex- 
antivortex pairbreaking effect. Experimentally it may 
manifest itself through a counterpart of Halperin-Nelson 
response which may be observable in systems with suf- 
ficiently strong dipolar coupling. In the second part of 
the paper we map DEC onto the tensionful vortex liq- 
uid state in layered superconductors. In that system we 
point out the emergence of a dipolar superfluidity which 
has a different origin because there is no interlayer pair- 
ing and carriers in layers have the same sign of electric 
charge. There the analogue of the dipolar superfluidity 
arises as a consequence of the fact that a molten lattice 
of composite vortices makes currents in individual layers 
dissipative while the counter-currents in different layers 
remain dissipationless. In these systems there is no need 
for interlayer pairing and carrier density is higher so the 
effective dipolar coupling may be relatively large. Besides 
that, the dipolar response can be used as an experimental 
tool to study vortex liquids in LSC, e.g. to unequivocally 
distinguish a tensionfuU vortex liquid state. 

We thank J.M. Speight for discussions and S. 
Shevchenko for communications and for informing us 
about a perturbative study of the weak dipolar coupling 
regime [^] and a study of Josephson coupling in DEC . 



This generalization of DEC may be called "multiflavor 
dipolar condensate" . 

In the both cases of DEC and LSC the dipolar super- 
fluidity will be destroyed by interlayer tunneling, but for 
small interlayer tunneling some of its physical manifesta- 
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